Dispersion of the neutron resonance in cuprate superconductors 
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We argue that recently measured downward dispersion of the neutron resonance peak in cuprate 
superconductors is naturally explained if the resonance is viewed as a spin-1 collective mode in a 
d-wave superconductor. The reduction of the resonant frequency away from the antiferromagnetic 
wave vector is a direct consequence of the momentum dependence of the d-wave superconducting 
gap. When the magnetic correlation length becomes large, the dispersion should become magnon- 
like, i.e., curve upwards from (7r,7r). 



PACS numbers: 71.10.Ca, 74.20.Fg, 74.25.-q 

Superconductivity and antifcrromagnctism arc two 
major phases of the high-T c cuprates. When viewed sep- 
arately, the corresponding ground states appear to be 
quite conventional: the parent compounds of high T c 
cuprates (such as La2Cu04 M) are exemplary Heisen- 
berg antiferromagnets, while overdoped cuprates in many 
respects resemble BCS-type d-wave superconductors [^|. 
However, full understanding of the interplay between the 
two phenomena is sorely lacking. Although most experts 
agree that antiferromagnetism is the ultimate cause of 
high-T c superconductivity, the intermediate steps are not 
yet clear ||. 

Perhaps the strongest experimental indication of the 
interplay between antiferromagnetism and superconduc- 
tivity in cuprates is the discovery of strong inelastic neu- 
tron scattering deep in the superconducting (SC) phase 
of materials with the highest T c . The most intense scat- 
tering at T « T c has been detected in YBCO |§ and 
Bi2212 H near the antiferromagnetic (AF) wave vector 
q = Q = (tt, it) and energy f2Q rs 40 meV. Experiments 
conducted with polarized neutrons indicate that the 
resonant scattering is due to electron spins. 

It is tempting to interpret the resonance as a magnon 
in a disordered Neel state, which becomes more visible in 
a superconductor. However, a magnon frequency 
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where c is a spin-wave velocity, clearly increases away 
from the AF wave vector Q = (Vr, tt) 0. Meanwhile, 
recent experiments on YBCO ||,[| have revealed the op- 
posite trend: the resonance energy decreases away from 
(tt, 7t) [see Fig. 0(a)]. In a disordered antiferromagnet, 
this can only be the case if spin response in the normal 
state is incommensurate ]To| . However, the data suggest 
H that, unlike in 214 materials pT] , the normal-state 
spin response in YBCO is commensurate. 

On second thought, the best superconductors of the 
cuprate family differ significantly from their parent AF 
compounds. They have a large, Luttinger-like Fermi sur- 
face in the normal state, and fermionic quasiparticles 



in the SC state, as shown by angle-resolved photoemis- 
sion spectroscopy (ARPES) |jl^]. Several authors have 
demonstrated that in this situation, an attractive ex- 
change force in the d-wave superconducting state binds 
an electron and a hole into a pair with total spin 1 [|l3] [l6| . 
This bound state is seen as a resonance in a spin response. 
The gross features of the neutron resonance, such as a 
non-monotonic variation of CIq with doping and persis- 
tence of the resonance in the pseudogap phase can be 
understood if the resonance peak is viewed as a collec- 
tive mode H-flT 
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FIG. 1. (a) A sketch of the resonance dispersion: experi- 
mental (solid line), a massive magnon in the disordered Neel 
state (dashed line), (b) Approximate shape of the Fermi 
surface in Bi2212. Filled dots are nodal points, where the 
fermions are gap less. The continuum of two-particle excita- 
tions at the corresponding wave vector qo = (n + 5, n + 5) 
(dashed arrow) starts at zero frequency. For lattice momen- 
tum Q = (n, 7r) (solid arrow) the continuum starts at a finite 
energy, (c) spin susceptibility at the wave vector (7r,7r). 

Formally, a collective mode has the dispersion given by 
Eq. (|l|). Its frequency is proportional to the static part 



of the inverse spin susceptibility (£ 
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However, this time the velocity c is not constant but de- 
pends rather strongly on the wave vector q. This de- 
pendence can be understood by noting that the bound 
state necessarily resides below the two-particle contin- 
uum [Fig. 0(c)]. For q along a zone diagonal, the con- 
tinuum starts at 2£ m i n (q) = |A(k)| + |A(k + q)|, where 
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both momenta k and k + q are at the Fermi surface, 
Fig. |l|(b). As q moves away from (tt, 7t), k and k + q 
shift along the Fermi surface towards the nodal points 
and the bound state inside the gap is pushed to lower 
energies. For q = q connecting the nodal points, the 
two-particle spectrum is gapless, i? m in(qo) = 0, so that 
the energy of the resonance must vanish (along with its 
strength). Obviously then, f2 q should decrease away from 
(tt, 7r), at least near qo. 

The above qualitative argument is quite robust and 
is based on two premises only: (i) the Fermi surface 
contains hot spots (i.e., points connected by the AF 
wave vector Q) and (ii) the superconducting gap has 
the d x 2_ y 2 symmetry. A more quantitative description 
of the resonance requires further assumptions. Most cal- 
culations of a spin response in a d-wave superconductor 
start with a model of free fermions on a square lattice, 
add an exchange interaction and employ the random- 
phase approximation (RPA) to compute spin suscepti- 
bility |l3|,[l4|,[l6| . The downward dispersion of the res- 
onance can then be obtained under favorable circum- 
stances Unfortunately, such calculations must be 
done numerically providing somewhat limited insight. 
Here we present an alternative, analytical approach to 
the problem. It is based on the observation that the 
behavior of a collective mode is determined largely by 
low-energy fermion degrees of freedom and is therefore 
insensitive to the physics at high energies. 

Our point of departure is a macroscopic spin-fermion 
model that describes low-energy fermionic quasiparticles 
(with a Fermi surface inferred from ARPES) interacting 
with collective spin fluctuations. The model is described 
by the effective action |l7j 

S = - fdr f 'rfT'^4 ff (r)G - 1 (k ) (T-r'))c k ,(r') 

+ \ f Q dr jf" dr' £ xo^q, r) S q (r) ■ S_ q (r') 

+5 r^5> q (r)-S_ q (T). (2) 
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where G *(k, r) = — Vk- (k— kp) is the bare fermionic 
propagator, and \o 1 ( < 1' r ) i s the bare spin susceptibility. 
This model is viewed as a low-energy version of a lattice, 
Hubbard-type model, obtained by integrating out high- 
energy degrees of freedom. Accordingly, it has a cut-off 
A < Ep and the bare xo (q> r ) is determined by fermions 
with energies exceeding A. If, as we assume, nothing 
special happens at high frequencies, Xo(q> the Fourier 
transform of XoiQ, t ) will have an Ornstein-Zernike form: 

Xo(q, fi) = Xo/(r 2 + (q - Q) 2 - 

Previous studies of Eq. (^|) focused on renormaliza- 
tion of the fermionic dispersion by spin-fermion inter- 
action, and on the form of the full x(Q,f2) [0. Here 



we consider the dynamical susceptibility at q ^ Q. To 
avoid unnecessary complications, we restrict calculations 
to momenta along the zone diagonal, q = (q, q). We also 
neglect strong coupling effects, assuming for simplicity 
that superconductivity is described by a d-wave version 
of the BCS theory. We have verified that strong-coupling 
effects (which modify fermionic propagator at low ener- 
gies) result in quantitative, but not qualitative changes. 

The full dynamical susceptibility %(q, ^) differs from 
Xo(q, fi) due to a bosonic self-energy n(q, il): 
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The static part of n(q, f2) (the contribution of the low- 
energy fermions to the inverse correlation length l/£) is 
small and, in fact, vanishes for linearized fermion disper- 
sion. On the other hand, the frequency-dependent part 
of n(q, f2) is substantially nonzero in the normal state. 
This is related to the fact that, for a Fermi surface with 
hot spots (as in YBCO and Bi2212) a low-frequency spin 
excitation with momentum q w Q can decay into two 
fermions at the Fermi surface (Fig. |l|). This gives rise to 
a universal relaxational term in n(q, Q): 



n(q, fi) w i Imn(q, fi) = i|fi|7q. 
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At small frequencies, n(q, f2) cx |fi| is much larger than 
the bare term 0(f2 2 ) in the susceptibility. It thus fully 
determines spin dynamics at low energies. The pref- 
actor 7 q is obtained by calculating the imaginary part 
of the particle-hole bubble. In two dimensions, 7 q = 
2g/(irv x v v ), where g = g 2 Xo is the effective spin-fermion 
coupling, and v = (v x ,v y ) is a (q-dependent) Fermi ve- 
locity at the hot spot kh s (q), whose components are de- 
fined as e k =v x k x + v y k y , e k+q = -v x k x + v y k y , where 
k = k — khs(q) Ordinarily, the q dependence is 

weak and can be neglected. However, for the wave vec- 
tor qo connecting the nodes, the velocities at kh s and 
khs + qo are antiparallel, therefore v y — and 7 q di- 
verges. We have verified that n(qo, f2) is of order yfi- 
Eqs. (ft) and (B) imply that in the normal state 
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where x — Q/j£ 2 , and q = q — Q. We see that spin 
response is (a) incoherent: no sharp peak in x"(q, ^) as 
a function of frequency f2; (b) commensurate: %"(q, J7) 
is peaked at q = Q at a fixed O. Both results are in 
agreement with the data . 

In a d-wave superconducting state, the bosonic self- 
energy is modified thanks to the opening of a supercon- 
ducting gap. It is now given by 
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Here A q is the fermion gap at one of the hot spots con- 
nected by q, i.e., A q = A(k), such that both k and k + q 
are at the Fermi surface. By virtue of the d x i_ y 2 symme- 
try, A(k) = — A(k + q), a condition we used in deriving 
Eq. (§). 

In the presence of a superconducting gap with the d- 
wave symmetry, Imn(q, f2) vanishes discontinuously for 

< 2|A q |. By virtue of the Kramers-Kronig relation, 
this discontinuity generates a nonzero Ren(q, fi) that is 
quadratic in f2 at low frequencies: 



n(q, Q) ~ 7 q 2 /A q as fl -> 0. 
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It is also essential that, for any q, n(q, 0) = 0: the open- 
ing of the d-wave gap does not change the magnetic corre- 
lation length. This result is not entirely surprising: spin- 
mediated d-wave pairing involves only fermions from op- 
posite sublattices and thus does not affect the correlation 
of spin within the same sublattice. 
Substitution of Eq. (0) into (||) yields 
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We see that in a d— wave superconductor the low en- 
ergy spin excitations are propagating, gapped magnon- 
like modes with the dispersion 
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Eqs. @-(||) are meaningful only if c q £ _1 < |A q |. Oth- 
erwise the use of a quadratic form for n(q, Q) is not jus- 
tified. Strictly speaking, near q = Q, c q £ _1 < A q only 
at sufficiently strong coupling, when fermion self-energy 
cannot be neglected. On the other hand, a quadratic 
frequency dependence of Ren(q, fl) at low fi is merely 
a consequence of a vanishing Imn(q, il) for fl < 2|A q |, 
where at strong coupling A q should be understood as 
the measured gap (i.e., a frequency where the spectral 
function has a 5— functional peak). We have explicitly 
verified that inclusion of strong coupling corrections into 
(0) only changes the overall factor in (Q). 

To verify that our analytical approach [based on a lin- 
earized fermion dispersion around kh s (q)] captures the 
essential features of the spin susceptibility, we present 
in Fig. U numerical results for the bare particle-hole 
susceptibility xo(q, ^) and its RPA counterpart x(q,fl) 
x(q, £1) = Xo(q. — ^Xo(q, ^)]- Both are calculated 

for a d x 2_ y 2 BCS superconductor with a tight-binding 
dispersion ek inferred from ARPES data |fj| . We expect 
that, at low enough n(q, f2) matches Xo(q, modulo 
a regular function of q. Similarly, x(q, fl) should agree, 
at low frequencies, with x(Qi 

In Fig. ||, we see that Imx (q, — for Q < 2|A q | 
and has a finite jump at 2A q , while Rexo(q, f2) is 
quadratic in ft at low frequencies, in agreement with what 
we obtained for n(q, fi). The peaks in Re%o(q, at 2A q 



are logarithmic singularities associated with the disconti- 
nuity of Imxo(2A q ). They are softened because we used 
complex frequencies with a small imaginary part. These 
singularities are indeed present in n(q, tt). A strong peak 
in the RPA susceptibility x(q, ^) is the resonance located 
at 1 = JRexo(q, This expression corresponds to 
n(q, fi) = + (q — Q) 2 in our analytical approach. We 
have also verified that the static susceptibility Xo^O) 
does not change between normal and superconducting 
states (i.e., £ is not renormalized). High-frequency fea- 
tures in Fig |^ — not not captured by our long-wavelength 
approximation — are lattice effects sensitive to the details 
of the fermion dispersion ek far from the Fermi surface. 
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2. Left: real and imaginary parts of the numeri- 
cally evaluated particle-hole bubble in a d— wave supercon- 
ductor, Xo(q, <^), at various q = (q,q), using a tight-binding 
fermion dispersion inferred from ARPES. The low-frequency 
features of xo(q, are in agreement with our analytical re- 
sults. Right: imaginary part of the corresponding full dy- 
namical susceptibility, Imx(f2, q). Observe that the resonance 
moves to lower frequencies when q 7^ Q = (tv, tv). 
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We now analyze Eq. (^|). If A q and 7 q were indepen- 
dent of momentum, c q would be constant and f2 q would 
be a converntional magnon-like dispersion (Q). In this 
situation, the spin resonance would remain commensu- 
rate and exist only at Q > J7q = cq/£. However, in a 
d-wave superconductor, A q decreases and 7 q increases 
when q deviates from Q. This effect accounts for the 
downward shift of f2 q . Furthermore, as c q vanishes at 
q^ q , the resonance frequency O q — > regardless of 
the spin correlation length £. Quite generally, near qo, 
Imx(q, should have a peak at a frequency O q < Qq. 

The behavior near q = Q is more complicated and de- 
pends on the values of £ and qo . We can roughly approx- 
imate the momentum dependence of c q as 



4 ~ 1 - (q/qo) 2 , 
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While the specific functional form we adopt here for c q is 
unimportant for capturing the qualitative behavior of the 
neutron resonance dispersion, the form adopted in Eq. 
( |Io| ) is appropriate for a magnetically mediated d-wave 
superconductivity in which A q is largest at q = Q |l7|] . 
Upon substituting this form into Eq. (m), we find 
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When qot; is small, the dispersion is negative for all mo- 
menta. Also, the residue of the peak in Eq. (||) scales 
aS Cq , 1 . e., it decreases and eventually vanishes as q ap- 
proaches qo- 

In Fig. H we plot f2 q ([ll]) and the intensity of the peak 
for go£ = 2. This value is consistent with optimally doped 
Bi2212 where |q | » 0.37r/a @) while £ = a - 2a jL9|. 
il q in ( pi] ) is rather flat near Q and rapidly drops away 
from Q. The residue of the peak at the downturn is 
already much smaller than its value at Q. Both these 
features are consistent with the data [Bl. 
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FIG. 3. Location and intensity of the resonance peak ob- 
tained from Fig. [] (left figure) and from Eqs. (|8]) and ( |ll| ) 
(right figure). 

For comparison, we also present in Fig. [| the disper- 
sion and intensity of the collective mode obtained in the 
RPA calculation, Fig. |[ It exhibits qualitatively similar 
behavior which we interpret as evidence that the disper- 
sion of the peak is insensitive to the details of the fermion 
dispersion at high energies. Note also that our estimate 
of the wave vector at which the energy and intensity of 
the neutron peak vanishes, qo « (0.87T, 0.87r), is roughly 
consistent with the corresponding value found in neutron 
scattering, q = (w, 0.8tt). Recall that, to simplify the 
analysis, we have only considered momenta along zone 
diagonal. 

Finally, we comment on the limit go£ ~^ 00 relevant to 
the underdoped side. In this limit, f2 q first increases with 
q and drops only very near q = qo . This implies that the 
resonance continuously evolves into the magnon of the 
AF state. It would be of interest to verify experimentally 
whether such evolution takes place. 

To summarize, we have demonstrated that the exper- 
imentally observed downturn of the resonant frequency 
H away from (it, 7r), accompanied by a rapid decrease of 
the peak intensity, occurs rather naturally when the res- 
onance peak is interpreted as a collective spin excitation 
in a d-wave superconductor. The unusual dispersion of 
the peak is related to a variation of the superconduct- 
ing gap along the Fermi surface. We have argued that 



the resonant frequency vanishes at a certain q = q that 
connects nodal points at the Fermi surface. Close to the 
commensurate wave vector Q = (jr, n), the dispersion de- 
pends on the spin correlation length: If the latter is small, 
|qo — QIC ^ 1) the resonance frequency f2 q decreases with 
|q — Q|. Conversely, when |q — Q|£ ^> 1, O q increases 
away from Q, reminiscent of a magnon dispersion in a 
disordered antiferromagnet. 
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